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ABSTRACT 

The heating of the ion-neutral (or ambipolar) diffusion may affect the thermal 
phases of the molecular clouds. We present an investigation on the effect of 
this heating mechanism in the thermal instability of the molecular clouds. A 
weakly ionized one dimensional slab geometry, which is allowed for self-gravity 
and ambipolar diffusion, is chosen to study its thermal phases. We use the 
thermodynamic evolution of the slab to obtain the regions where slab cloud 
becomes thermally unstable. We investigate this evolution using the model of 
ambipolar diffusion with two-fluid smoothed particle hydrodynamics, as outlined 
by Hosking & Whitworth. Firstly, some parts of the technique are improved to 
test the pioneer works on behavior of the ambipolar diffusion in an isothermal 
self-gravitating slab. Afterwards, the improved two- fluid technique is used for 
thermal evolution of the slab. The results show that the thermal instability may 
persist inhomogeneities with a large density contrast at the intermediate parts 
of the cloud. We suggest that this feature may be responsible for the planet 
formation in the intermediate regions of a collapsing molecular cloud and/or 
may also be relevant to the formation of star forming dense cores in the clumps. 



Subject headings: ISM: clouds - Hydrodynamics - ISM: magnetic fields - diffusion 
- methods: numerical - ISM: evolution. 
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1. Introduction 

Observations of molecular clouds show the existence of small density fluctuations threaded 
by magnetic fields. The recent progress of the molecular cloud observations has established 
that the small and tiny scale structures are very ubiquitous. For example, Langer et al. 
(1995) observed small condensations with size from 0.007 to 0.021pc in Taurus Molecular 
Cloud 1. The mass of these fragments is estimated to be < 0.01— 0. ISM©. Studies of the time 
variabihty of absorption hues indicates the presence of fluctuations on scales of 10~^pc (5-50 
AU) and masses of IQ-^Mq (Boisse et al. 2005). At larger scales (about 10,000 AU), Pan 
et al. (2001) find significant differences in CN, CH, and CH+ absorption hues. The SPITZER 
experiment has begun producing higher spatial resolution mid-infrared maps (Churchwell et 
al. 2004) and revealed the fine structure of the star forming regions. The advent of major 
new facilities in the coming years should yield several breakthroughs in this field (Andre 
et al. 2008). A higher resolution observation in future is assuredly expected to reveal hid- 
den, small structures, as well as sub-stellar objects with very small mass (Atacama Large 
MiUimeter Array, ALMA, becoming partly available in 2011, fully operational at 2013). 

The classical problems of cloud fragmentation are relevant to star formation because they 
show that there is a preferred scale (Jeans length) of gravitational instability as soon as one 
considers an initial structured region with some length-scale (see, e.g., Larson 1985). By the 
way, it is not easy for the small and tiny scale structure to form as a result of gravitational 
instability, since their size are much smaller than the Jeans length. According to Langer 
et al. (1995) indeed, the small scale structures appear to be gravitationally unbound. This 
suggests that some fragmentation mechanisms other than pure Jeans gravitational instability 
may be important in the clouds. Molecular clouds are thought to be supersonically turbulent, 
since they exhibit supersonic linewidths (Zuckerman & Palmer 1974). Base on the idea of 
turbulence, some mechanisms to generate substructures inside molecular clouds have been 
proposed (Elmegreen 2007). As well, a further process based on magnetohydrodynamic 
waves have been suggested to produce inhomogeneities in the clouds (Fohni et al. 2004). 
Another basic physical process that could trigger formation of density fluctuations is thermal 
instability, because it acts on time-scales that can be much shorter than the duration of the 
turbulent motions (Gilden 1984). The effect of thermal instability in fragmentation of the 
clouds is clearly a competitive process relative to other mechanisms, although all could 
operate. 

In molecular clouds, a neutral molecular gas is intermixed with an ionized component 
that is tied directly to the magnetic field. Birk (2000) has used the two-fiuid technique to 
find the thermal condensation modes in weakly ionized hydrogen plasma. Nejad-Asghar & 
Ghanbari (2003) studied the effect of linear thermal instability in a weakly ionized magnetic 
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molecular cloud within the one-fluid description. Nejad-Asghar (2007) has recently made the 
assumption that the molecular cloud is initially an uniform ensemble which then fragments 
due to thermal instability. He find that ambipolar drift heating is inversely proportional to 
density and its value, in outer parts of the cloud, can be significantly larger than the average 
heating rates of cosmic rays and turbulent motions. His results show that the isobaric 
thermal instability can occur in intermediate regions of the cloud; therefore it may produce 
the cloud fragmentation and formation of the condensations. 

The study under the one-fiuid approach did not precisely consider the effect of the ion- 
neutral friction, since they did not treat the flow as two fluids which arc composed of ions 
and neutrals. Then, Fukue & Kamaya (2007) revisited the effect of the ion-neutral friction 
of the two-fluid on the growth of the hnear thermal instabihty. Their results indicate that 
the friction with the magnetic field affects the morphology and evolution of the interstellar 
matter. Falle et al. (2006) described numerical calculations using an AMR magnetohydrody- 
namic code that show that thermal instability may have an important role to play in the 
formation of the hierarchical structure of molecular clouds. Here, we use the two- fluid nu- 
merical simulation to study the evolution of weakly ionized molecular clouds in the nonlinear 
thermal instability regime. 

Many authors have developed computer codes that attempt to model ion-neutral dif- 
fusion. Black & Scott (1982) used a two-dimensional, deformable-grid algorithm to follow 
the collapse of isothermal, non-rotating magnetized cloud. The three-dimensional work of 
MacLow et al. (1995) treats the two-fluid model in a version of the ZEUS magnctohydro- 
dynamic code. An algorithm capable of using the smoothed particle hydrodynamics (SPH) 
to implement the ambipolar diffusion in a fully three-dimensional, self-gravitating system 
was developed by Hosking & Whitworth (2004, hereafter HW). They described the SPH 
implementation of two-fluid technique that was tested by modeling the evolution of a dense 
core, which is initially thermally supercritical but magnetically subcritical. 

In this paper, we flrstly re-formulate and improve the two-fluid SPH implementation of 
the ambipolar diffusion in an isothermal self-gravitating one dimensional slab. Afterwards, 
we use this two-ffuid SPH technique to investigate the nonlinear thermal evolution of the 
slab. In this way, the continuum equations of the self-gravitating slab and its isothermal 
case are given is section 2. The two-ffuid SPH technique is presented in section 3. Section 4 
devoted to the chosen physical scales, initial conditions, and the computer experiments to 
explore the temperature proflle and density fluctuations in the contracting one dimensional 
molecular cloud. Finally, section 5 is allocated to the summary and conclusions with some 
prospects. 
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2. Gas dynamics 

In the absence of reliable information about the structure of the molecular clouds, it 
is reasonable to examine simple geometries such as sphere, disks, cylinder, etc., in order to 
gain some insight into the evolution process. Here, we consider an one-dimensional lightly 
ionized molecular gas with purely transverse magnetic field, thus, all variables are functions 
of distance z to the center and time t only. 



2.1. The continuum equations 

The exact fraction of the total fluid that is ionized depends upon many factors (e.g. the 
neutral density, the cosmic ray ionization rate, how efficiently ionized metals are depleted 
on to dust grains). Here, we use the expression employed by Fiedler & Mouschovias (1992), 
which states that for 10* < n < lO^^m"^, 

P. = 6(//^ + eV'), (1) 

where in standard ionized equilibrium state, e ~ 7.5 x 10~^^kg^/^.m~^/^ and e' ~ 4 x 
10~^^kg^/^.m~^^''^ are valid. In reality, the gas in this case is very weakly ionized, thus, 
we adopt the approximation p = p„ + pj ^ p„ in fluid equations. In this case, the molecular 
cloud is considered as global neutral which consists of a mixture of atomic and molecular 
hydrogen (with mass fraction X), helium (with mass fraction Y\ and traces of CO and other 
rare molecules, thus, the mean molecular weight is given by 1/p = X/2 + y/4. 

We write the continuity equation of neutral part as its common form 

while, the relation ([1]) is used to determine the ion density wherever it is required in the 
two-fluid equations. The momentum equation of neutral part then becomes 

dv 1 9 , 52 

at p oz 2po 

where the gravitational acceleration g obeys the poisson's equation 

I = -4vrGp, (4) 
and the pressure is given by the ideal gas equation of state 

P = —pT (5) 
p 
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where R is the molar gas constant and the temperature T is approximated the same as for 
both neutral and ion fluids (Tj — Tn — T). 

The calculations are usually further simplifled by the assumption of an isothermal equa- 
tion of state. This last assumption is, however, unnecessarily crude since molecular gas is 
expected to cool and heat rapidly and appropriate cooling and heating function may be 
estimated (e.g.. Goldsmith 2001). The thermal energy per unit mass is generally given by 

where the mean internal (rotation and vibration) energy of an H2 molecule is included. At 
temperature below ~ 200K, the rotational and vibration degrees of freedom of molecular 
hydrogen are not excited, so it actually behaves like a monatomic gas. The energy equation 
follows from the flrst law of thermodynamics, that is 

du pdv 

where ^{p,T) is the net cooling function 

^^(p,T) = ^w(y^)^^"^ - i^CR + Tad), (8) 

where Tcr and Tad are the heating rates due to cosmic rays and ambipolar diffusion, re- 
spectively, and A(ji) and are the parameters for the gas cooling function that here we 
use the polynomial fltting functions, outhned by Nejad-Asghar (2007), as follows 

( tW -0 --8-98-0-87aog-) 
VJ.kg ^s V ^0 

-0.14(log-)^ (9) 
riQ 

= 3.07 - 0.11(log — ) - 0.13(log —Y, (10) 

no no 

where no — lO-'^^m"^. The heating of cosmic rays is given by (e.g.. Goldsmith 2001), 

FcR ~ 3.12 X 10-^ J.kg-\s-\ (11) 

while the heating due to ambipolar diffusion can be physical derived by considering the drag 
force (per unit volume) as follows 

Vad = (12) 
P 

where 

fd = lADPipyd: (13) 
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where 7ad ~ 3.5 x 10^''m^.kg~^.s~^ represents the colhsion drag coefficient. 

The magnetic fields are directly evolved by charged fiuid component, as follows: 

where the last term outlines the ambipolar diffusion effect with drift velocity, 

lAD(^PPiOZ 2/io 

which is obtained by assumption that the pressure and gravitational force on the charged 
fluid component are negligible compared to the Lorentz force because of the low ionization 
fraction. 



2.2. Isothermal molecular layer 

In isothermal case, the thermal energy per unit mass (E]) is overplus, and the pressure 
in momentum equation ([3]) then becomes p = a^p where a is the isothermal sound speed. 
Following the many previous treatments, a further simpliflcation is possible if we introduce 
the surface density between mid-plane and z > as 



p{z',t)dz'. (16) 



By transformation from {z,t) to (a, t), the drift velocity is given by 



1 ^(^\ 



and the equation (iMj) becomes 

d Id, dB. 



dt p ■jAD^Poda p^/'^il + e' p'^/"^) da 
With the above, fleld equation (j4]) can be integrated to give 

g = -AnGa, (19) 

while the equation of continuity ([2]) and the equation of motion ([3]) take the form 



dz _ 1 
da p 



(20) 
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and 

|_._4.G.-^(aV+|;). (21) 

respectively. The slab is assumed to be in quasi-magnetoliydrostatic equilibrium at all times, 
supported against its own self-gravity by the magnetic and gas pressures. The loss of flux 
from ambipolar diffusion is exactly compensated for by the compression of the slab which 
is necessary to maintain equilibrium. In this approximation, the left-hand side of equation 
fl2T|) is zero, and we may integrate the force balance to obtain 

^ + aV = 27rG(a^-a^) (22) 

where integration constant cToo is the value of a at z = oo (where p is zero). 

Following the work of Shu (1983), we introduce the non-dimension quantities 
(T _ ^ 27rG(Too f; B 



(Too' SttGo-^ ' ' 2aooV7i"/ioG 

~ 2V27rG' 27rG'o-oo. , . , . 7ADe ^^^n 

7ADe a (27rG)5/V^ aV27rG 

so that we rewrite the basic equations (ITSl) . (I2U]) and fl2^ as follows: 

-(^) = -^( gg) (24) 



52 + p=l-a^ (25) 
a£ _ 1 

and the drift velocity (fT7|) as 



(26) 



1 ^5 

= + (27) 

Here, a natural family of initial states is generated by assuming that the initial ratio of 
magnetic to gas pressure is everywhere a constant, i-e., -B^/p = ao at t = 0. Then one 
finds from equations (!25l) . (126|) and (p7|) that 

P(2,i=o) = 727^7 — 7' (28) 

cosh (2:/ Zoo) 



2ao av27rG smh(z/zoo) 
VI + "0 7ADe 1 + e'po cosh^(2/zoo) 
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where po = 27rGcr^/a^(l + ao) is the central density of the slab at t = 0, and 2:00 = 
a a/ (1 + ao)/27rGpQ is a length-scale parameter. 

Figure [T] shows the initial neural density, ion density, and drift velocity in the cloud and 
outercloud medium for po = 3.8 x lO^^^kg.m"^, a = 0.55km.s^^, and = 1- Shu (1983) 
solved the equations (p^ - (p71) by finite difference techniques under the initial and boundary 
conditions 

B = I "° I (1 - ^2)1/2 at t = and 5 = at a = 1. (30) 
\1 + aoj 

The integrations are carried out to time t = 20 with initial conditions ao = 1 and 10, and 
the results are presented graphically in his paper (Shu 1983). As the magnetic field leaks 
from the neutral gas, the volume density of the neutrals shifts in profile from equation fl2Sl) 
to that case with ao = 0, and the drift velocity gradually settles. 



3. Numerical scheme 

Since the molecular gas is expected to cool and heat rapidly, considering the thermal 
energy is appropriate. For this aim, we produce a two-fluid SPH implementation more 
efficient than HW, to study the thermal evolution of the self-gravitating one dimensional 
molecular cloud. In the technique of HW, the initial SPH particles are represented by two 
sets of particles: magnetized ion SPH particles and non-magnetized neutral SPH particles. 
For each SPH particle we must create two separate neighbor lists: one for neighbors of the 
same species and another for those of different species. Consequently, each particle must 
have two different smoothing lengths. In the following sections we refer to neutral particles 
as a and (3, and ion particles as a and b; the subscripts 1 and 2 refer to both ions and neutral 
particles. 

The neutral density in place of neutral particles is estimated via usual summation over 
neighboring neutral particles 

Pn,a = ^rnf3Waf3, (31) 

while in place of ions, pn,a, is given by interpolation technique from the values of nearest 
neighbors. The ion density is evaluated via equation ([1]) for both places of ions and neutral 
particles. In the particle approximation of SPH method, the infinitesimal volume in the 
integrations at the location of particle a is replaced by the finite volume of that particle 
AV^, which is related to the mass and density as = paAVa- Since the ion density is 
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evaluated by equation ([T]), no by the usual summation rule 



Pi,a = '^mbWab, (32) 

/3 

we must update the mass of ion a as follows 

= (33) 

ra 

in each time step so that the above summation/interpolation routine for ions ( l32l) being 
accurate. 

The SPH form of the drift velocity of ion particle a is given by HW as 
1 . 1 „ dWab ^ dWa, 



Vd,a 



—BbBa-T— - 2^m„n„ft-— ]. (34) 



lADPn,a PoPi,a ^ Pi,b dZa ^ dZa 



where Uab is the usual artificial viscosity between ion particles a and b (Monaghan 1992). 
The artificial viscosity is reconsidered in SPH to prevent inter-particle penetration, unwanted 
heating, and unphysical solutions. Nejad-Asghar, Khesali Sz Soltani (2008) has recently 
considered the coefficients in the Monaghan's standard artificial viscosity as time variable, 
and a restriction on them is proposed such that avoiding the undesired effects in the subsonic 
regions. Here, we use the Monaghan's standard artificial viscosity, since the cloud contraction 
is quasi-hydrostatic and there is not supersonic motions and shock formation during this 
contraction. Keeping in mind the second golden rule of SPH which is to rewrite formulae 
with the density inside operators (Monaghan 1992), we can optimize the drift velocity of 
HW as follows 



Vd,a 



lADPn,a '2^PQPi,a ^ Pi,b " dz, 

nib T-r dWa 



— H^b— — , (35) 

^ Pi,b dZa 

where two extra density terms are introduced, one outside and one inside the summation 
sign. This comes as a result of the approximation to the volume integral needed to perform 
function interpolation. The drift velocity at neutral places are used to estimate the drag 
acceleration 

0'drag,a = lADPi,aVd,a, (36) 

instead the method of HW who used the expression of Monaghan & Kocharayan (1995). 
Since there is no any analytical expression that allows us to calculate the value of drift 
velocity in place of the neutral particles, we use the interpolation technique that starts at 
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the nearest neighbor, then add a sequence of decreasing corrections, as information from 
other neighbors is incorporated (e.g.. Press et al. 1992). 

In the usual symmetric form, the self-gravitating SPH acceleration equation for neutral 
particle a is 

= 9c.~2^T^p{-^ + -o+^o.p)-T— + adrag,o. (37) 

where Qa is the gravitational acceleration of particle a. The ion momentum equation assum- 
ing instantaneous velocity update so that we have 



Va = y2 —'"P^ap + Vd,a (38) 

where the first term on the right-hand side gives the neutral velocity field at the ion particle 
a, calculated using a standard SPH approximation. 

The SPH equivalent of the energy equation ([7]) is 

dUa 1 Pa Vl3 ^ „ . dWap ^ 

^ = 9 2^ ^M-^ + -2 + n„/3)t;„/3^r— - ^a. (39) 
at Z ^ pp UZa 

The temperatures of neutral particles are calculated from equation ([6]), and the temperature 
of any ion is assumed to be same as its nearest neutral neighbor. Finally, the magnetic 
induction equation in SPH form is replaced by 

dt ^ Pb dZa 

where the usual notations of the ion fluid are used. 



4. The computer experiments 

The chosen physical scales for length and time are [/] = 200AU, and [t] = lO^yr, 
respectively, so that velocity unit is approximately [v] = Ikm.s"^. The Newtonian constant 
of gravitation is set G = l[m]~^[/]^[t]~^ for which the calculated mass unit is [m] = 4.5 x 
lO^^kg. Consequently, the derived physical scale for density, energy per unit mass, and drag 
coefficient are [p] = 1.7 x lO'^^kg.m'^, [u] = lO'^J.kg-^ and 'Jad = 1-8 x 10^^[lf[m]-^[t]-\ 
respectively. In this manner, the numerical values of e and e' are 1.8 x 10^^[/]^^/^[m]^/^ and 
3.5 X 10^^'^[/]^^^/^[m]^/^, respectively. The magnetic field is scaled in units such that the 
constant po is unity. Since the magnetic fiux density has dimensions 

" [t][chargey ^^^^ 
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while /io has dimensions 

specifying /ig = 1 therefore scales the magnetic field equal to [B] = S.lnT. With aforemen- 
tioned units, the thermal energy per unit mass (El) is represented by 8.3 x 10-=^(5X/4+3y/8)T, 
the heating rates due to cosmic rays and ambipolar diffusion are Tcr = 1.1 x 10^'^[n]/[t] and 

TAD,a = lADPi,o,v\a^ (43) 

respectively, and the parameters for the gas cooling function are 
log = -4.48 - 0.87(log ■ 



u]l[t]J ■ ■ ' "^2.24 X 10-4' 



-0.14(log — T)^ (44) 

^ ^ 2.24 X 10-4^ ' ^ ' 



/?. = 3.o7-o.ii(iog ^^/" 



2.24 X 10- 
2.24 X 10-4' 



-0-13(log ,,/" ,_J l (45) 



4.1. Initial setting 

The initial conditions for this simulation are a parallel magnetic field directed perpen- 
dicular to the z-axis so that the initial ratio of magnetic to gas pressure is everywhere a 
constant (ao = 1), and a density profile given by the equation (1251) . The magnetic field is 
assumed to be frozen in the fluid of charged particles and the central density is assumed to 
be po = 2.24 x 10-4[p]. We choose a molecular cloud which has a mass fraction of molecu- 
lar hydrogen and helium X = 0.75 and Y = 0.25, respectively, and has an initial uniform 
temperature of Tq = 50K. We assume that the cloud slab is spread from z = — 78[/] to 
z = -1-78 [/] (according to Fig. [T]). The initial values of the cooling and heating functions are 
shown in Figure [2l As presented in this figure, the isobaric thermal instability criterion, 

OA dT 

is satisfied in the outer parts of the cloud, thus, these regions are thermally unstable while 
the inner part is stable as outlined by Nejad-Asghar (2007). 

Implementation of the boundary particles in the diffusion processes is an important 
problem in two-fluid SPH simulation. In ambipolar diffusion process, the ion particles are 
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physically diffused through the neutral ffuid, thus, the ions will be bared in the boundary 
regions of the cloud (i.e. without any neutral particles in their neighbors). Both the cloud 
and boundary regions contain ion and neutral particles. The complete system is represented 
by discrete but smoothed SPH particles (i.e. Lagrangian sample points) so that coagula- 
tion of particles and fragmentation of the cloud is truthfully revealed. Since it is desirable to 
have initially the same numerical resolution for both fluid components, we use N/2 ions and 
N/2 neutral particles. We set up boundary particles (4/ii up and down in z) using the hnear 
extrapolation approach (from the values of the inner particles) to attribute the appropriate 
drift velocity, drag acceleration, pressure acceleration, energy rate, and the magnetic induc- 
tion rate to the boundary particles. We check the position of ions before making a tree and 
nearest neighbor search, so that we do not consider the bared ion particles in the boundary 
regions of the simulation at next time-step. 

The present SPH code has the main features of the TreeSPH class so that the nearest 
neighbors searching are calculated by means of this procedure. The selection of time-step. 
At, is of great importance. There are several time-scales that can be defined locally in the 
system. For each particle 1, we calculate the smallest of these time-scales using its smallest 
smoothing length, hi, i.e. 

Ati = Cco^j^minf-^, (47) 

l^^ll VA.l Cs,l 

where va = B / ^JJI^pl is the Alfven speed of ion fluid and Ccour is the Courant number which 
in this paper is adopted equal to 0.3 (for numerical stability). The evolution were carried 
out to time < u > / < A, F >~ 0.43/0.025 = 17.2[t] so that the fragmentation of the cloud 
via thermal instability may be revealed. 



4.2. Results 

The molecular cloud has an initial uniform temperature of To = 50K. The Fig. [2] 
shows that the cooling is greater than heating at initial time t = 0, but by passing the 
time, the drift speed of ions in the outer layers increases, thus, the heating rate growths 
in those regions. Increasing of the heating rate in the outer layers causes to stabilize those 
regions, and the intermediate layers become thermally unstable. The evolution of cooling 
and heating rates, causes to increase the temperature of the slab at outer regions while 
decrease the temperature at the intermediate parts of it. The temperature profiles at times 
t = 3.5 [t], 10.5 [t] and 17.2 [t] are shown in Figure [31 This figure shows that the outer layers 
of the slab heat up, the intermediate layers cool down, and the inner layers hardly change. 
It is obvious that the instability of the cloud at the intermediate parts, causes the formation 
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of two relative cool regions in those areas. 

Since, in course of time, the thermal instability is shifted to the intermediate regions of 
the slab, the condensation mechanism is occurred there as depicted in Fig. m The growth 
of thermal instability results a density imbalance between the preformed condensations and 
their adjacent medium. As a result of this instability process, the particles are dragged 
from the adjoining of the condensations so that the density fluctuations growth and the 
fragmentation of the parent slab into cold, dense, low-mass cloudlets may be appeared. To 
compare the effect of the thermal instability, we firstly perform the isothermal simulation 
that have initially T = 50K. The thermal and density evolution of the slab are carried 
out to times t = 3.5 [t], 10.5 [t] and 17.2 [t]. The isothermal density profile and the density 
comparison with the isothermal case, are shown in Figure |H This figure shows the resulting 
isothermal density evolution at the top panel while the density contrast, relative to the 
isothermal case, is shown in the bottom panel. In this run, the thermal instability is fully 
developed and makes density fluctuations in the intermediate parts of the cloud. 



5. Summary and conclusions 

Molecular clouds have a hierarchical structure that extends from the scale of the cloud 
down to much smaller masses for unbound structures. A weakly ionized self-gravitating one 
dimensional slab geometry is assumed in order to gain some insight into the fragmentation 
process. The two-fluid SPH technique was used to investigate the nonlinear thermal evolution 
of the slab. The initial conditions for this simulation are an inverse square cosine-hyperbolic 
profile for initial density and a parallel magnetic field directed perpendicular to the slab. 

The initial values of the cooling and heating functions are shown in Figure [21 According 
to this figure, the isobaric thermal instability criterion is satisfied in the outer parts of the 
cloud. In course of time, the thermal instability is shifted to the intermediate regions of the 
slab so that the condensation mechanism is occurred there, thus, these regions are thermally 
unstable while the outer part is stable. The evolution were carried out to time 17.2 [t] so that 
the fragmentation of the cloud via thermal instability is revealed. 

The cooling and heating rates cause to increase the temperature of the clouds at outer 
regions while decrease the temperature at the inner part of it. It is obvious that the instability 
of the cloud at the intermediate parts, causes the formation of two relative cool regions in 
those areas, which is shown in Figure [31 The rapid growth of thermal instability results a 
density imbalance between the cloud and the surroundings. The isothermal density profile 
and density contrast, relative to the isothermal case, are shown in Figure [H In this run. 
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the thermal instability can fully develop and makes density fluctuations in the intermediate 
parts of the cloud. 

This feature may be responsible for the planet formation in the intermediate parts of a 
collapsing molecular cloud and/or may also be liable for the formation of star forming dense 
cores in the clumps. The macro-velocity fields in the molecular clouds are highly turbulent 
and supersonic. It is then of uppermost importance to consider both collision and merger 
of the formed condensations. Merging is possibly the main onset mechanism to form dense 
cores, which likely evolve to star formation. 
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Fig. 1. — The initial drift velocity, ion density, and neutral density in the cloud and outer- 
cloud medium (gray region) for po = 3.8 x 10~^^kg.m~^, a = 0.55km.s~^, and = 1- 
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Fig. 2. — The initial values of the cooling (solid) and heating (dash) functions versus position 
and neutral density. 
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Fig. 3. — The temperature profile versus position at times t = 3.5 [t], 10.5 [t] and 17.2 [t]. The 
initial temperature of the slab is chosen as uniform which is shown by dash line. 
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Fig. 4. — The isothermal density profile versus position at times t = 3.5 [t], 10.5 [t] and 17.2 [t] 
(top panel), and the density comparison with the isothermal case at those times (bottom 
panel). The dash line in the top panel is the initial central density of the cloud. 



